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COMBUSTION THEORY FOR LIQUIDS WITH A FREE SURFACE 
III. SPECIAL PROBLEMS 

S. N . Mil'kov, G. 3. Sukhov, L. P. Yarin 

The quasi-one-dimensional theory developed in [1,2] can be /3 * 
used to solve a broad group of problems related to combustion 
of liquids with a free surface. We will illustrate this with 
examples of approximate calculation of two processes which are 
interesting in terir.3 of technical application: extinguishing a 

flame when a liquid in a container burns and when a liquid from 
a porous layer burns. 

A flame above a free surface can be extinguished by rapid ! 

H 

mixing [3] = The mechanism of this process involves the y 

following. As the liquid is mixed, the free surface's 

temperature and, therefore, evaporation rate drop. If heat 

removal from the flame is nearly constant, the decrease in flow 

of vapor to the front reduces the reaction rate and, 

ultimately, results in extinction critical heat exchange i 

conditions [4] . j 

To find the values of parameters corresponding to 
extinction, we will use a one-dimensional model of a flame 
(fig. 1) and a quasi-heterogenous diagram of a flame front ; 

[5] . Let us write the equations for material and thermal 
balance at the front 

Gfr-o — G^+o = ( (1) 

qiG^o “ £4+0) * =s (?*+» <?*-<» ( 2 ) 

where and are, respectively, diffused flows of 

vapors approaching the flame front and departing from it; 


*Numbers in the margin indicate pagination of the foreign text. 
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q M and ty+o 

Xexp i-E/RT*) , 


thermal flows at the front; W 4 , -* spict* X 
the rate of effective quasiheterogenous reaction.* 



Fig. 1. Temperature 
distribution during 
combustion of a liquid 
with a free surface, 
la, b - Gas medium; 

2 - Reacting liquid; 

3 - Substrate. 


Flows incorporated into equations (1) and (2) can be found 
by solving the system of quasi-one-dimensional transfer 
equations followed by calculation, of the derivatives according 
to known temperature and concentration profiles [1] . To 
simplify tne computations, in future approximation equations 
for diffusion and thermal flows we will use 



is the nominal thickness of the layer of oxidizer above the 
free surface, at whose boundary c sc . 


Let us introduce the coefficient for vapor combustion 


*The remaining values in (1) and (20) and subsequent equations, 
as well as subscripts, are taken from [ 1 , 2 ]. 
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efficiency 


<P — 1 ~ v ' i 

where ctm-cwCl-- £*/&,) is vapor concentration at section £”&t> 
if there is no combustion. 

Substituting equations (3) and (4) into (1) with regard for 
(7) yields 


q>» 



( 8 ) 


x-z/Z>-U{b ~U)fo‘, e«**RTjE. 


Considering , we obtain . Assuming that 

parameter P«i2, during steady combustion until flame-out, we 
obtain an approximated equation for [1] . 




Where p = fic6o; fi=>Q/C(,*. Using this relationship and clapeyron- 
Clausius' equation p =*=Qxexp {— qjR/R b EQ 0 ), vib express T as an 
exponential function of the free surface's temperature 


t = t 0 [1 —x exp (— g r n /0o)] _1/s . 

» 


Here T o 


( 2 g ; -^) 1/2 ln(i + Q); q a = quR/ER b . 


Transforming (2) using (3) -(7), we obtain* 


<P 2 


(Q^-Qq) Si g 


(9) 


Since the scalar concentration of vapor c^ is 


*In deriving (9) , it was assumed that Le=l. 


sufficiently 
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close to c b0 , we can, assuming that 

C bm » X OXP (-Sa/0o)| 

give (9) the following form: 

<p s — oxp (? b /©oK 


( 10 ) 


( 11 ) 


whe re g =* qRn/cpiE.- 

To determine the free surface's temperature 0 Q , one must 
use the equation for heat balance at the phase interface. Heat 
flow incorporated into these equations reflects both heat 
release from gas to liquid and. transfer of heat by thermal 
conductivity and convection in the liquid. The solution to /5 

the problem includes calculation of transfer processes in the 
gas and liquid if there is forced convection (mixing) in the 
latter. However, to analyze the effect of mixing on the heat 
mode of combustion, there is no need to solve such a 
complicated problem. It can be limited merely to considering 
the relationship between combustion temperature (reaction rate) 
and free surface temperature, since physical considerations 
make it clear that mixing the liquid results in a reduction in 
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Pig. 3. Temperature distribution 
during combustion of a liquid from a 
porous layer. 

1 - Gas phase; 2 - Porous substance 
(lb*, region of vapor filtration 
through the porous layer) . 


Therefore, regarding as a parameter, let us consider 
the path of curves «p^(0^) and variables <f>-6 (fig. 

2). The function <j > 2 ( 6 ^) is described by an S-curve starting 
at origin and, as 0 ^, having an asymptote cp = 1 . Increasing 
(acclerating mixing in the gas flame) shifts the curve 
<^( 0 ^) to origin. The function of> 2 (0^ ) is linear,, A change 
in shifts it along the x-axis and abruptly changes its 
slope. This reflects two consequences of the change in free 
surface temperature: a change in heat removal from the front, 

which is proportional to G^-Gq, and vapor concentration at 
the phase interface, which is an exponential function of 0 ^. 
Figure 2 clearly shows that a reduction in 6 ^ causes 
combustion temperature to drop and vapor concentration at the 
flame front to rise. At a certain critical 0 Q , curves 
°Pj_(G^) andtppfGjj) lose their common point at the upper 
temperature level, which corresponds to flame-out (point H). 


Using ( 8 ) and (11) , we will find the critical free surface 
temperature. Assuming that curves (0^ ) andq^tQ^) are 
tangent, we derive an equation which relates critical 
combustion efficiency and critical temperature 

•P-p-l-e^re- (12) 

e 4> ^ e 6 

Equating the right sides of (11) and (12) , we find critical 
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temperature as a function of free surface temperature and 
adjusted calorific value q: 



( 1 ?) 


The minus sign before the radical in (13) represents ^6 

combustion; the plus sign, extinction. Prom (8) and (11), we 
have 


*(e„) • 

t(Q 0 ) + e m * 


4- (0* - 0 O ) exp (f?n/©o)’ 

9 • 


(14) 


Equality (14) is vaJ id for any steady state, including critical. 


Solving the system of equations (13) and (14) , one can find 
the critical free surface temperature 0g as a function of 
assigned parameters q and q p . 

Let us consider the problem of combustion of a liquid from 
the porous layer of a substance. Experiments show that two 
different liquid combustion modes can be employed, depending on 
the system's hydraulic characteristics [6], Given low porous 
layer hydraulic resistance (high medium porosity, low reacting 
liquid viscosity) , the filtration supply mechanism meets the 
requirements of the reaction. Equilibrium is established in 
the system, and the evaporation front is localized at the 
layer's surface. If the medium's permeability diminishes, the 
equilibrium is disrupted, and the evaporation front shifts to 
within the layer. Its outer surface heats to a temperature 
above the liquid's boiling point. Therefore, the liquid- and 
gas-phase flow field can be divided into four typical zones 
(fig. 3) . In the first, 2, the reactive liquid, exposed to 
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capillary forces and pressure applied from outside, is filtered 
through a layer of porous substance of thickness ^><$2* ftChe 
liquid warms due to heat supplied to the porous substance from 
the flame front, in region lb', vapor formed during 
evaporation is filtered through a layer of porous substance of 
thickness 4^^* In zones lb and la, the vapor mixes with the 
oxidizer found in the environment, and the components interact 
chemically. 

Xn a quasi-one-dimensional approximation, diffusion 
combustion of a liquid filtered through a layer of porous 
substance is written as follows in a system related to the 
boiling front* 





Pi u i if” 8 * ~' a<tC 

(15) 

for 

the 

region 

Ifil < i < +oo ( 





dc * r~\ d c 

Pi u i"df ~ 9lD d£ ’ 

(16) 




dT[ ’ rdX 
,9l u l c pl~tfr 5=1 

(17) 




u l dp'- 1 ' rpt 

“ pl ~ r 

(18) 

for 

the 

region 






dT„ iT z 

mpjUjiCpg — %2 jgt < 

(19) 





(20) 


for the region where m and k^ are porosity and the 

*At Le=l, one diffusion equation is sufficient to describe 
heat- and mass-transf er processes in the gas above the porous 
layer . 
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fl,ltca ( fcion factor? and ^ 2 ' Q ^^ QC ^ VQ coefficients of 
thermal conductivity for the filled porous medium in zones lb 1 
and 2. 

Boundary conditions which are taken into account during £]_ 

integration of system (15) -(20) take the form 


£->- + 00 : 

T-+T wt 

c — >■ 0, 


T m T$, 

c — c at> 

S ■" ?«• * 

T-T* 

P *■* PlVr, 

6-0: 

T ** 

~ C m Cq( 

£ das! 


P ~ Pi*- 


( 21 ) 


Balance equations reflecting the preservation of heat and 
substance flows toward typical surfaces should be added to 
conditions (21) : 


i-U, h (-3r) 1+I - - ®i D [-if)*-,'' 

(22) 

i-u *.(^L,-4SrL’ 

(23) 

. ’ ,, fdr'A iar,\ 

6-°, 

(24) 

- Pl w l c M Pl& ("jf) “* Pl^l — Pl u l m t 


f * Q \ 

(25) 

|i-Q*oxp( 

(26) 


However, without fully solving the problem, one can identify 
certain essential features of this process on the basis Of the 
simplest possible model.* For example, it is of practical 


*Integrating system (15) -(20) does not involve any particular 

difficulties. 
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interest to evaluate limit pressure drop A/> Bp »• p 2 * -- p 1+ , 
at which the evaporation front shifts inward into the porous 
layer. I'o find it, assuming _ m 2V»2V,. 2»"*2V, 

' 1 HS 


(xL+o^ (%L. 




1 


'2V, 

1 K 


61 


' (5 )-s ^ > 


(h is the depth of the porous layer; t v9 , the liquid's 
boiling point) , we obtain a system of three algebraic equations 


\ JV. _ A ? 




1 HS 


A’ r 6i 


■T, 


KS 


»Wh 


■), 


Hi 


®6l 

' T 




2* 


k 0 b P 




ki/h *+ 


h ~~ §61 
Pi 


(27) 

(28) 
(29) 


which define T , ^ , and u 2 as a function of thermal, 
filtration, and geometric characteristics of the subject 
system. Specifically, solving (27) -(29) relative to Ap, we find 


Ap 


S&1 ~ Pi) Pi ' 1 ( \ 

.WiV* rif * A ”6#1 /’ 


(30) 


introducing ^^“0 into (30* , we obtain the desired equation, 
which relates limit values for parameters representing the 
transition of combustion from a mode in which the liquid 
evaporates from the porous material's surface to a mode where 
the surface of the phase transition lies within the porous 
substance layer 

(31 > 


Assuming that the layer of porous substance is sufficiently 
thick and the thermal effect of the chemical reaction is high 
(W* «0/£2 t a T U2 tt'T 2 +, we obtain the following 

Ap ap « aji/mhc^t, • p u /£2 • O/On, * (32) 


where pu^Pi/p*; a, — l|/c P1 p,. Equation (32) shows that, all other 
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conditions being equal, the limit gradient is governed by the 
ratio of reaction heat to vapor-formation heat. 

Using data irom [6], we can approximate /f*A/>, p , which 
corresponds to the movement of the boiling front to the outer 
surface of the porous layer. For example, for kerosene in 
sodium chloride (h«120 mm, d*16 mm, m*0.4), the transition from 
filte .ed combustion (boiling front within the layer) to free 
combustion (evaporation from the layer's surface) occurs when 
the grain cross section changes from 10 to 340 microns. The 
values k^Ap*0.1S *10'’ 2 and 2.16*10“^, which correspond 
to these conditions, also determine the range of possible limit 
values km&puv , Calculation with (32) and data from (3) over 
the length or the flame yields k*&p BV =0.09 *10 -5 m 2 /sec, 
which is close to the lower limit of this range. 
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